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An exactly soluble non-linear interaction Hamiltonian is proposed to study fundamental 
properties of the entanglement dynamics for a coupled non-linear oscillators. The time-evolved 
state is obtained analytically for initial products of two coherent and two number states and 
relevant informations are extracted from the dynamics of various quantities like subsystem linear 
and Von Neumann entropies, quadrature mean values, variances and Q-functions. We determined 
the re-coherence time scales and found among the interaction terms present in the Hamiltonian 
the one responsible for the entanglement in both cases. We identify the existence of two regimens 
for the entanglement dynamics in the case of initially coherent states: the short time, phase 
spread regimen where the entropy rises monotonically and the self-interference regimen where the 
entropy oscillates and re-coherence phenomenon can be observed. We also found that the break 
time from the first regimen to the second one becomes longer, as well as the re-coherence and re- 
versibility times, as the Planck's constant becomes much smaller than a typical action in phase space. 
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I. INTRODUCTION 



Nowadays it is widely accepted that quantum entangle- 
ment is an essential ingredient for the implementation of 
quantum information processing devices . It is notable 
that until recently quantum information science was re- 
stricted over discrete, finite dimensional Hilbert space 
elements, however, interest on continuous variables has 
been developing in various contexts: from teleportation 
[H, cryptography Q to cloning Q. This is closely con- 
nected to the recent advances of enhancing nonlinear cou- 
pling via electromagnetically induced transparency (KIT) 
mechanism using a Bose-Einstein condensate, which 
has opened possibilities of strong non-linear interaction 
of ultra slow light pulses [|| of tiny energies belonging to 
different modes of electromagnetic fields ||^, ^ . Such ad- 
vances ranges from building a quantum logic gates based 
on photonic qubits without resorting to photonic crys- 
tals ^ , entangling continuous variable states ^ , and 
preparing entangled states of radiation from mixed ther- 
mal states 1 11 in conventional media. Facing these recent 
developments, it is important to understand the entan- 
glement properties of Kerr-type interaction particularly 
for Gaussian states, which are one of the states that can 
be generated by simple experimental devices like beam 
splitters and phase shifters. 

In the present work, we have studied the entanglement 
properties of a model describing a bipartite system of two 
degrees of freedom, an integrable version of two quartic 
oscillators ||l2|, |l3j, coupled via a Kerr-type interaction. 
One could think of some coupled cavities scheme with a 
Kerr medium |l^ as a possible realization of these 
type of interactions. The nonlinearity of each of the sub- 
system oscillator keep some of the properties present in 
one degree of freedom such as collapses and revivals of 



the quadrature mean values of each field. Also, this sol- 
uble model in which the nonlinearity is present can be 
used to understand the role of Kerr-type nonlinearity on 
the entanglement dynamics of initially coherent Gaus- 
sian states. Another case we have considered is the non- 
classical number states for the initial disentangled 
state, which does not entangle via Kerr-type coupling, 
but via usual bilinear coupling in a rotating wave ap- 
proximation (RWA), and show several differences in the 
entanglement dynamics. 

The paper is organized as follows: In section II the 
model is introduced and the exact solutions for initially 
disentangled states, in the cases of number states and co- 
herent states, are presented. Analytical expressions for 
the subsystem density operators and some mean values 
are calculated. Section III is reserved to present the exact 
subsystem linear entropies for both initial states, analyze 
the conditions for re-coherences and recurrences and dis- 
cuss the differences in the entanglement dynamics in the 
proposed cases. Also, we examine the differences of the 
short time dynamics where no interference phenomenon 
is possible and the longer time dynamics where such ef- 
fects do affect the evolution of the subsystems. Another 
issue treated here is the limit where Planck's constant 
becomes much smaller than the typical action in phase 
space. As we shall see, the two cases present distinct 
semi-classical behaviors. Finally in section IV we present 
our conclusions. 



II. THE MODEL 

Our theoretical model will be inspired on two field 
modes in a cavity with a low-loss Kerr media, described 
by the creation and annihilation operators and d\ 
{k — 1,2), such that each mode has the usual non-linear 
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interaction term of the form X^^k^l: fill- Besides this 
non-Unear self-interaction term we add two couphng in- 
teraction between the field modes: one is the usual RWA 
coupling [|l7[ a nd the other is the Kerr-type non-linear 
interaction |15| , 
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The resulting Hamiltonian can be written in terms of the 
number operators Nk = ^^.A^, of the new oscillators as 
follows: 



N1+N2 + I 



(5) 



Using Eq.(U), we can connect the two different basis of 
the two oscillators Hilbert space by the relation 



where the constants g and g pq | depend upon the third- 
order non-linear susceptibility x ^^'^ losses will be 
considered. Since our purpose here is to have an exactly 
soluble model with the nonlinearities in such a way that 
we could trace out the effect of each interaction term on 
the entanglement dynamics, we will just set g = 2g and 
also choose the resonant case {oji — u!2) such that the 
Hamiltonian can be re-written in the following form: 
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Notice that the entire Kerr type interaction is now in 
the form of the square of the free Hamiltonian Hq, and 
hence the following commutation relations hold 



0. 



(2) 



This property is important for two reasons: first, it 
is associated to the presence of a constant of motion 



{J\f = a[ai + 0,20,2)', and second, it allows us to study 
separately the action of each interaction term in the evo- 
lution operator on the initial state, and hence its conse- 
quences on the entanglement process. In what follows, 
we study the cases of initially uncorrelated number and 
coherent states. 



A. Analytical solution for product of number states 

For the case of initially disentangled number states of the 
harmonic oscillators: 



(3) 



Wc can write H in Eq.(^) in terms of two new quar- 
tic oscillators diagonalizing the RWA-coupling, using the 
transformation proposed by Zoubi et al [ |l9| , 
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where 
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(7) 



The sub-index 'a' or 'A^ indicates the bosonic represen- 
tation the ket belongs to. It is to be noticed that at the 
right hand side of Eq.(^, only those states with the total 
number Af fixed by the initial state are present. Using 
the Hamiltonian in the diagonal form, Eq.(|5[), and Eq.(|^) 
the time-evolved state can be written as 



(8) 



where the quantity $ = uJot{U+l) + XW+hgt{J\f+ 1)^ 
is a glo bal p hase. This phase factor will be relevant in 
Section II B, where this result will be used to construct 



the solution for the case of coherent states. Hence, for 
the number states, Eq.(H) shows that the dynamics of 
the stationary states are completely determined by the 
RWA coupling and the effect of non-linear part of the 
Hamiltonian is only in the global phase. 

After a little algebra we can re-write the above state 
Eq.(0) in a more instructive form: 



\m)^r^u\t)f'^^i'\t) |o,o), 



fi(")/j.\_ (a* cos At-ia[, sin At)" 



(9) 



The sub-indexes "12" ("21") of the operators f empha- 
size the entanglement features of the dynamics. Here, 
we omitted the remaining global phase, and used the 
fact that the vacuum state in both representations spaces 
must be the same. 

Now, we can extract some informations of entangle- 
ment process of number states under the action solely of 
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the RWA interaction. At times T; = j{l + I), with I 
integer, Eq.(|9|) gives 



_j)ni+n2(_i)i(«l+"2) |„2,ni) , (10) 



indicating that the system is in a disentangled state, but 
not the same as the initial one. However, at times t; = 
I J we recover the exact initial state (modulo an overall 
phase) : 



IV'(tO) = (-i)'("^+"^)K,«2) 



(11) 



Hence, for the case of initially disentangled number 
states, we have found two characteristic times: (a) the 
re-coherence times Ti, for the which the subsystems re- 
cover the purity of the initial state, and (b) the recur- 
rence times Ti, when the evolved state becomes equal 
to the initial state. This allows us to classify li^iO)) in 
Eq.(H) as a reversible state for this particular interaction. 
Consequently, we may define a period of reversibility for 
such states: th = ^. For the special case of equal ini- 
tial number states (ni =712), the difference between the 
re-coherence and recurrence times disappears and the re- 
versibility occurs earlier, as can be seen in Eq. ([l0|) . 



B. Analytical solution for coherent states 

Consider an initially disentangled product of two co- 
herent states as follows: 
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The time-evolved state can be obtained using the previ- 
ous result Eq.(0), including the A/'-dependent phase $ 
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where a-k{t) = afee"*''^''^^''^)*. In this result we have 
already omitted a global phase. 

It is easy to see that for g = the RWA-coupling does 
not entangle the oscillators, since the two summations of 
Eq.(|l3|) can be factorized, and each oscillator remains in 
a coherent form 
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where 

(ii{t) = (aie" 



^'^«*cosAi-za2e-''^''*sinA0, (15) 



and bk[Pk{t)] = e(''''°'^ ''^°'') is the displacement opera- 
tor in the phase space of the fc-th oscillator. 



Using this particular result, we can solve the general 
case of initially disentangled coherent states in a more 
intuitive way. The commutation relations (||) allow us 
to apply separately the piece of the evolution operator 
associated to the non-linear term of Eq. |l|, and use 
the previously derived results [Eqs.( |l^, ^^]. The final 
exact solution for the temporal evolution of two initially 
coherent states is given by 

|f/)(t)) = e~^e ^ \ai,a2) 

= e^|/3i(i),/32(i)) 
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Here, the expressions of /3fc(i) are the ones depicted in 
Eq. (|l|). 

In order to calculate other quantities like the mean 
values and variances of the quadrature operators, and 
discuss the phenomenon of collapses and revivals, let us 
calculate the density operator of the system. Since we are 
interested in the case where the global system is isolated, 
the total density operator is a projector onto the state 

\m) 
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where we have defined A = -t- |/32| , a constant of 
motion associated to the mean value of J\f . We shall 
denote from here on (3k{t) simply as (3k- We calculated 
the reduced density operators pk{t) (A: = 1, 2) by tracing 
over the undesired degree of freedom, corresponding to 
one of the original oscillators: 
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with a similar expression for p2 (t). The field quadrature 
operators Qk and Pk are given in terms of the creation 
and annihilation operators as follows: 
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and analytical expressions for the mean values are given 
by 



Qk){t) = Trk[QkPk{t)] 



(20) 



= v^Rek(t)e-3'"«*e-A(i-^"'"''')l , 
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'Pk){t) = Trk[PkPk{t)] 



(21) 



Pk{t)e 



where we defined the frequency associated to the non- 
linear term as LOg = hg. The above expressions contain 
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FIG. 1: (a) Collapse and revival phenomenon in the evolution 
of the expectation value of the quadrature Qk and (b) the 
respective variance. Analytical results for ^ = 20, ^ = 2, 

Qfeo = Pko ~ 1.0, A — 4.0 and h — 1. The same qualitative 
behaviors are observed for {Pk)(t) and (APfe)(t). 



low integer fractions of the re-coherence times. But, in 
contrast to the case coupled to a reservoir we will 
see in what follows that at the revival time the subsys- 
tem entropy goes to zero (re-coherence), and the system 
do come back to the initial state (recurrence) under ap- 
propriate conditions. Another initially disentangled state 
that can be solved is the product of the form \ai) (8) |n2), 
where this type of analysis can be done and find the times 
at re-coherence allows superposition states [p2|. 



III. ENTANGLEMENT PROPERTIES AND ITS 
SEMI-CLASSICAL BEHAVIOR. 

In this section we are going to discuss the entan- 
glement dynamics of time-evolved states found in Sec- 
tion H; namely, the product of number states and co- 
herent states. We analyze the subsystem entropy, one 
important tool that gives a measure of entanglement for 
globally pure bipartite systems, and also the Husimi dis- 
tribution or Q-function, in order to follow the evolution 
of the partial distribution in phase space. The Q-function 
for the global system is usually defined by: 



Qiq^p)^ - ill 72I p(i) I71 72) 

TT 



(22) 



where we choose — "^y^*" (fc = 1,2). Particularly, we 
investigate the role played by the non-linear interaction 
term in the Hamiltonian (Q) which appears associated 
to the frequency Ug in the various analytical expressions 
derived. Moreover, we would like to find the behavior of 
the subsystem linear entropies in the semi-classical limit. 



Entanglement properties for number state: 
linear and Von Neumann entropies. 



exponentials of oscillatory terms, which produces the col- 
lapses and revivals phenomena at times odd multi- 
ple of 2^ and integer multiples of ^ respectively, as 
shown in Fig.|l|(a) where we plotted the mean value of 
Qk quadrature as a function of ujgt. 

One can also calculate analytically the variances 
(A(5/;)(i) and {APf^){t). The results are complicated ex- 
pressions which we shall omit. However, we illustrate 
its general behavior in Fig.0(b), where we can see the 
same periodic structure present in the quadrature mean 
value temporal behavior. These results for the subsys- 
tem variances are very similar to those obtained for the 
one degree of freedom quartic oscillator |12 . The most 
instructive information is the relation between the col- 
lapse and the phase space's delocalization in phase of the 
state, as can be seen in the projected Q-function that will 
be illustrated in the next section for a particular set of 
parameters. 

Due to the fact that here we have two interacting sub- 
systems, their entanglement prevent the appearance of 
superpositions of coherent states for each oscillator for 



We are interested in the calculation of both, the sub- 
system linear entropy (SLE) and the Von Neumann en- 
tropy, usually defined by 



4 = 1 - Trkipl) 
Sk = -Trkipk'^npk), 



(23) 



where the label k is associated to one of the degrees of 
freedom. For the initially disentangled number states 
Eq.(||), the reduced density operator is diagonal in the 
number state basis; therefore, it can be directly expressed 
in terms of its eigenvalues as follows 



ni+n2 
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(t)^ ^ Mt)\i){i\, 
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where the trace condition ^ A;(t) = 1 must be satisfied. 
The eigenvalues Xi{t) are given by: 
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(26) 



Then, in terms of the reduced density operator eigen- 
values, the subsystem entropies assume a famihar form: 



ni+n2 



Skit) 



6k{t) = i~ J2 ^'W' 

- J2 A,(t)ln(A,(<)). 



(27) 



The dependence of the eigenvalues on the periodic func- 
tions confirms the rule for the re-coherence times, Eq.(p^' 
and Eq. ( [TTI ) . To illustrate what has been discussed above 
we present some simple examples: 



• For 1-0(0)) = |1,0) we have: 

[ Sk{t) = - cos^ Xtln (cos^ \t) - sin^ At In (sin^ \t) 

• For |-0(O)) = we have: 
Skit) = 2142^(5 + 3 cos 4At), 

Sk{t) = - cos^ 2At In (cos^ 2\t) - sin^ 2At In 



(28) 



(29) 



• For 1-0(0)) = |2,0) we have: 
Sk{t) = 2i£ip(i3 + 3cos4At), 

•S'fe it) = - cos" Xt In (cos** \t) — sin" Xt In (sin* Xt) + (30) 
_£i£i2Min| 



These analytical results for the SLE are shown in 
Fig.^(a) . An interesting fact is the appearance of the os- 
cillatory structures between two successive re-coherences 
for particular values of Uk ■ This is illustrated in Fig.||(b) 
where the SLE is plotted for ni — 3 and several suc- 
cessive values of n2- We can see that both, the number 
of oscillations and the maximum value of SLE, increases 
as 712 is increased. This feature can be understood from 
Eq. (p4|) , where the number of accessible states is given 
hy ni + n2 + I which increases with increasing 712. This 
accounts for the increase in the maximum value, and also 
the subsystem is allowed to pass through many more pos- 
sible mixed states. Notice that no dependence on h or 
g is essentially left in the expression of the time-evolved 
state (H), and the same happens to the calculated SLE's 



[Eqs.(|28[§^)]. This is the case where the non-linear inter- 
action plays no role, and the semi-classical limit is related 
to the presence of many photons (large Uk)- 
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FIG. 2: (a) Subsystem linear entropy for three different sets 
of initial number states: |1,0) (solid line), |1, 1) (dashed line) 
and 1 2,0) (solid thick line), (b) Subsystem linear entropy for 
initial number states (|ni , 712)), with 712 = 3 and several values 
of Til (legends). Note the particular behavior for 711 — 772 at 
At=|. 



Entanglement properties for coherent states 
and semi-classical behavior. 



In order to study the entanglement dynamics and re- 
versibility properties for the initial product of coherent 
states, we calculate the SLE defined in Eq.( |2^). Us- 
ing previous result Eq.(p^, an exact expression for this 
quantity can be calculated: 

S,it) = 1-Tniplit)) 

Z ^ -rn I 777 I 



X e 



' [hgt(n-~m)] 



(31) 



In a similar fashion as in the case of the number states, 
we can find the re-coherence times studying the condi- 
tions under which Siit) is equal to zero, indicating that 
the subsystems are disentangled. A simple analysis of 
Eq.(^ll) shows that this happens in the following situa- 
tions: 

1. For any initial conditions excluding the vacuum 
state (|0, 0)), the re-coherence times associated to 
the non-linear interaction are hgTi = In, when the 
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argument of the sine- function is an integer multiple 
of tt: 



T; = ; — integer). 
hg 



(32) 



2. For those initial conditions such that one of the 
arguments /3fc(i) = 0, one can find other instants 
when SLE is zero. These conditions are associated 
exclusively with RWA coupling and special initial 
conditions. In terms of the quadratures of the 
initial values a^, using the expression ( |l5| ) one can 
write these conditions as follows: 

9192 + P1P2 = 0, 
and we have disentangled states for times p^ : 

h = j-arctan{-^), or ta = jarctan{^). 



Among these re-coherence times, we can identify those 
at which we also have recurrences by examining Eq.(16). 
Only the first class of re-coherences produces recurrences, 
and this happens when the following condition is satis- 
fied: for a given value ofTi, when the ratio and l-^ 
are integer numbers. 

Now we will choose a particular parameter sets to illus- 
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FIG. 

A 



3: Subsystem linear entropy 5k as function of loqI for 
0.2, ^ = 0.1, qko = PkQ = 1.0, A = 4 and several values 
of 7^ = 1: 7^ = J (dotted line), 7?- = | (dotted thick line), 
TZ = (dashed line), TZ = ^ (solid line) and TZ = ^ (solid 
thick line). 



trate the semi-classical limit of the SLE, since in contrast 
to the case of the number states, we have an explicit de- 
pendence on h, actually on the frequency LUg = hg. This 
allows us to study the semi-classical behavior of the SLE. 
The natural parameter to measure the 'quantumness' of 
the syst em h ere is the ratio TZ = j^, where A (defined in 



Section [IB) is a characteristic action in phase space. In 
Fig.^ we plotted the time evolution of the SLE for several 
values of TZ. Here, we adopted the convention to fix the 



value of A and vary h instead of the opposite way. The 
first thing to be noticed is the fact that all curves coin- 
cides at the short time scale, where the SLE increases 
monotonically until it reaches the maximum value. The 
maximum value of SLE depends on TZ, which increases as 
we let 7?. ^ 1. This has to do with the increasing num- 
ber of accessible states as we let the spectrum become 
denser, and also implies in loss of information. We will 
call this first regimen the ''phase spread regimen' in con- 
nection to the behavior of the Q-function of each subsys- 
tem that we shall see in what follows. After reaching the 
maximum value, oscillations start to happen in the SLE, 
which can be seen as a partial recovering of coherence, 
until the first re-coherence time Ti = ^ (see Eq.(|2|)) at 
which (5fe = and the subsystem recovers purity. This 
second regimen will be called ^self-interference regimen' 
since the time-evolved subsystem Q-function shows the 
phenomenon of self-interference typical of the Kerr-type 
nonlinearity. This is consistent with the fact that in the 
limit TZ ^ Q the re-coherence time goes to infinity and 
the initial purity will never be recovered. We will call 
break time ti,, the time at which the transition from the 
phase spread regimen (rising) to the self-interference (os- 
cillating) regimen occurs. It is clear from Fig.|^ that this 
time increases with the inverse power of TZ, and we also 
expect it to go to infinity in the classical limit . 
Let us illustrate the differences between the two regimens 
by means of the subsystem Q-function which will show 
the proper signatures in each regimen. In Fig. ^ we plot- 
ted Q-function in various instants of the phase spread 
regimen for the same Hamiltonian parameters used in 
Fig. ^, and the quantumness parameter TZ — 0.025 . 
The sequence of plots is in the quadrature plane of the 
oscillator-1, beginning at < = until t k, tb- Where the 
center of the initially coherent Gaussian wave-packet fol- 
lows essentially the classical trajectory, as predicted by 
the Ehrenfest theorem , circulating around the origin 
and, due to the nature of the self-interaction term, the 
packet itself spreads in phase angle in the phase space. 
During the interval of time before the front of the packet 



reaches its tail |^ 
In Fig.||, anot 



no re-coherences can happen, 
ler sequence of contour plots of Q- 
functions of the subsystem- 1 shows a time evolution dur- 
ing the self-interference regimen in the interval of time 
after the break time th, until the first re-coherence time 
Ti . We also add a plot at the recurrence time (ri = 2Ti 
for this particular set of parameters). It is remarkable 
the appearance of several peaks along the annular re- 
gion, a signature of self-interference phenomenon where 
a kind of standing waves with M-peaks forms at times 
t — Ti/M (for t > tb), when the SLE assumes the value 
« 1 — -i. at a local minimum. This kind of behavior is 
a hallmark of the self-interference regimen, an essentially 
quantum phenomena which is at the core of re-coherence 
and reversibility in this system. This sequence of phe- 
nomena phase spread and self -interference is reproduced 
many times latter since the entanglement process in this 
case is reversible. We call the attention to the fact that 
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IV. 



CONCLUSIONS 
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FIG. 4: Contour plot of subsystem Q-function on plane 
{qi,pi) during the phase spread regimen, for initial conditions 
Qio = Pio = q2o = P20 = 1-0, tJo = 1, A = 0.2, g = 0.1 
and TZ = 0.025. The first plot correspond to the initial time 
(upper-left corner) and the subsequent times are chosen in 
the interval < t < tt- 
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FIG. 5: Time evolution of contour plots of a subsystem Q- 
function during the self-interference regimen for TZ = 0.025 in 
the plane (gi,pi), with the same parameter values of Fig. ^ 
(a) t = Ti/8, {h)t = Ti/7, (c) t = Ti/6, (d) t = Ti/5, (e) t = 
TlR/i, (f)f = Ti/S, (g) f = ri/2, (h) t = Ti (re-coherence) , 
(i) t — 2Ti — Ti (recurrence). 



in this particular case shown, the standing waves are not 
Schrodinger cat states, in fact it is a mixture of cat states, 
in contrast to the case reported in Rcf. ^ for a sim- 
ilar Kerr-type Hamiltonian. We present this calculation 
in the Appendix |A| . 



We have solved analytically the problem of two reso- 
nant RWA-interacting fields in the presence of non-linear 
Kerr-like interactions. The time-evolved state was ex- 
actly determined and this allowed us to identify analyti- 
cally how the collapses and revivals are produced in the 
quadrature mean values for the initially coherent state. 
All properties of the entanglement dynamics have been 
studied for initial states which are products of both co- 
herent and number states. 

In particular, we presented all the necessary conditions 
to the re-coherence of the initially non-entangled number 
and coherent states. We have calculated the exact ex- 
pressions for the subsystem entropies in both cases. Also, 
some conditions for the recurrence has been established 
for the coherent state case, related to the commensura- 
bility of the physical frequencies of the model. 

We also identified, in the case of coherent initial states, 
two distinct regimens of entanglement: the first one 
(phase spread regimen) happens during the time where 
the initial coherent state spreads in phase angle in the 
phase space; whereas the second one (self-interference 
regimen) occurs when the phase spread state starts to 
self-interfere. The time at which self-interference be- 
comes important for the evolution of this type of ini- 
tial states, we call it break time (tb), and it delimits the 
beginning of the essentially quantum processes responsi- 
ble for the re-coherence. We also have shown that the 
self-interference of each oscillator produces the standing 
waves, where no single Schrodinger cat state is allowed 
but a mixture of Schrodinger cat states, consistent with 
the entanglement. Finally, we show how in the semi- 
classical limit (h 0) both the re-coherence times and 
the recurrence times go away and the entanglement pro- 
cess becomes irreversible for all practical purposes in this 
type of model for Gaussian initial states. 
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APPENDIX A: A DEMONSTRATION OF 
MIXTURE OF SCHRODINGER CAT STATES 



We first re- write Eq. (16) in a more compact notation 
m)) = ^c„(7i)c„(72)e-'^''*("+'")'|n,m)(Al) 
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c„(7) = e-^^ (A2) 

where 7^ = /3fce~*^^''* and (3k are functions of time. We 
are interested in the instants tr « = - = Ti - , where r 
and s are mutually prime with r < s. Using the discrete 
Fourier transform 

9=0 

where 

^ r s, if r is odd and s is even, or vice- versa, 
1 2s, if both r and s are odd; 

and 

a('-'^) = i^e-"K'=^-2f), (A5) 

fc=0 

we re- write Eq.( |Al[ ) at the mentioned instants as follows: 
l^(^r,.)> = ^c„(7i)c„(72)e-"5("+'")' |n,m) 

n.m 

= ^ QgUp ^ c„ (771, 6"*^™ = ) c„i (?72p) |n, m) 

g,p^O n,m 

(A6) 



= ^ UqUp ^ |?7i, e"*^'""^ ) (g) c„i (?72p) |m) 

l-l 

= ^ |C„) (g) ^OpCm (?72p) |m) . (A7) 

m p=0 

Here rjkq = 7fce^'^'^T = /?fce~*^'^(^^~) , and we have 
defined the following cat-like states for the oscillator- 1 

i-i 

\Crn) = 

Then, constructing the global density operator and trac- 
ing over the oscillator-2, i.e., summing over the left over 
number states, we finally get the following mixture of cat 
states: 

Pl{tr,s) = ^U|C™)(C„| (A9) 

m 

l-l 

Cm = X! ^p"'*p'Crn{V2p)c*^{V2p')- (AlO) 
p,p'=0 
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